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ABSTRACT 


The South Galactic Cap u-band Sky Survey (SCUSS) is a deep u-band imaging survey in the Southern Galactic 
Cap, using the 90Prime wide-field imager on the 2.3 Bok telescope at Kitt Peak. The survey observations started in 
2010 and ended in 2013. The final survey area is about 5000 deg” with a median 5o point source limiting 
magnitude of ~23.2. This paper describes the survey data reduction process, which includes basic imaging 
processing, astrometric and photometric calibrations, image stacking, and photometric measurements. Survey 
photometry is performed on objects detected both on SCUSS u-band images and in the SDSS database. Automatic, 
aperture, point-spread function (PSF), and model magnitudes are measured on stacked images. Co-added aperture, 
PSF, and model magnitudes are derived from measurements on single-epoch images. We also present comparisons 
of the SCUSS photometric catalog with those of the SDSS and Canada—France—Hawaii Telescope Legacy surveys. 


Key words: catalogs — surveys — techniques: image processing — techniques: photometric 


1. INTRODUCTION 


The South Galactic Cap u-band Sky Survey (SCUSS; see X. 
Zhou et al. 2015, in preparation for more detail) is a deep 
imaging survey in the South Galactic Cap in u-band with an 
effective wavelength of 3538A. It is an international 
cooperative project between the National Astronomical Obser- 
vatories of China and the Steward Observatory of the 
University of Arizona. The survey utilizes the 2.3m Bok 
telescope at Kitt Peak. The camera, installed at the prime focus, 
provides a field of view (FOV) of about 1 deg”. The adopted 
filter is similar to the u band of the Sloan Digital Sky Survey 
(SDSS; York et al. 2000). The SCUSS project started in the 
summer of 2009 and began its observations in fall 2010. 
Survey observations were completed in the fall of 2013. The 
final survey area is about 5000 deg? with uniform imaging 
depth, of which more than 75% is covered by the SDSS 
footprint. The median imaging depth for point sources is about 
23.2 mag at a signal-to-noise ratio (S/N) of 5 with a 5 minute 
exposure time. Table 1 gives the basic characteristics of 
SCUSS. 

The main goal of the survey is to supply input photometric 
catalogs to select spectroscopic targets for the Large Sky Area 
Multi-Object Fiber Spectroscopy Telescope (Cui et al. 2012). 
In addition, by combining with other bands in large-scale 
photometric surveys, such as the SDSS and Panoramic Survey 
Telescope & Rapid Response System (Pan-STARRS; Kaiser 
2004), the survey data can be used for a wide range of scientific 
investigations, such as Galactic structure, Galactic extinction, 
galaxy photometric redshift, galaxy star formation rates, and 
stellar populations of nearby galaxies. 

This paper describes the data reduction pipeline specially 
designed for the SCUSS survey. There are some instrumental 
issues that need to be specially handled, such as issues in 
overscan, substructures in bias, and crosstalk. In addition to 
detecting sources on SCUSS images, we also provide 
photometry for SDSS objects using consistent object para- 
meters. Section 2 introduces the SCUSS survey and related 


facilities. Section 3 describes the basic image processing. 
Astrometric and photometric calibrations are presented in 
Sections 4 and 5, respectively. Section 6 provides image- 
quality statistics of the observations. Section 7 describes image 
stacking. Sections 8 and 9 present photometry methods and a 
comparison with other surveys, respectively. A summary is 
given in Section 10. 


2. THE SURVEY AND FACILITIES 


SCUSS is a u-band imaging survey in the northern part of 
the southern Galactic Cap. The survey originally covered the 
region of Galactic latitude b < —30° and equatorial latitude 
6 > —10°, with a total area of about 3700 deg”. It was further 
extended to the Galactic anti-center region and the extra area 
covered by the SDSS, with a final survey area of about 
5000 deg? (X. Zhou et al. 2015, in preparation). Normally, 
there are two exposures for each field and the total exposure 
time is 5 minutes, which generates images 1-1.5 mag deeper 
than the SDSS u band. In the following, we give a general 
description of the telescope, camera, and filter that were used in 
SCUSS. 


2.1. Telescope 


The Bok telescope” is a 90 inch (2.3 m) telescope operated 
by the Steward Observatory of the University of Arizona. It is 
located on Kitt Peak, whose latitude is +30°57/46”5 and 
longitude is 111°36’01”6W. The elevation is about 2071 m and 
the typical seeing is about 1”5. The telescope runs year-round 
except on Christmas Eve and during a maintenance period in 
August. The accuracy of the absolute pointing of the telescope 
is recorded as 3” over the entire sky. 


5 http: //james.as.arizona.edu /~psmith/90inch /90inch.html 
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Table 1 
Survey Summary 


Telescope 2.3 m Bok telescope 

Site Kitt Peak in Arizona 

CCD 2 x 24k x 4k CCD array 
Field of view 1°08 x 1°03 

Filter u (3538 A) 

Integral time 300 s 

Magnitude limit 23.2 mag 

Survey area ~5000 deg” 

Observation period 2010-2013 


2.2. Camera 


An imaging system, named 90Prime, is deployed at the 
prime focus (corrected focal ratio: f/2.98; corrected focal 
length: 6829.2 mm). The detector is a CCD array consisting of 
four 4k x 4k backside-illuminated CCDs. They are STA2900 
CCDs made by Semiconductor Technology Associates, Inc. 
and backside processed at the University of Arizona Imaging 
Technology Laboratory. These CCDs have been optimized for 
the u-band response, giving a quantum efficiency at u band of 
about 80%. Figure 1 displays the layout of the CCDs on the 
focal plane. The edge-to-edge FOV is about 1°08 x 1°03. The 
pixel scale is 0”454. There are inter-CCD spacings in the 
center of the array: 166” in right ascension and 54” in 
declination. Each CCD is read out by 4 amplifiers located at the 
corners. Each detector has 4096 x 4032 physical pixels and 20- 
row pixels of overscan for each amplifier. The full well is about 
90,000 electrons or 65,000 data numbers (DNs). The current 
system gain is set to be about 1.5 electrons per DN. The dark 
current is about 7 electrons per pixel per hr. The readout time is 
about 30 s and the average readout noise is about 8.8 electrons. 
The response non-uniformity with the u filter, which is the 
standard deviation divided by the mean of a u band flat field 
image, is about 2%. Table 2 summarizes these parameters of 
the detector. In 2010, CCD #4 had problems and only a 
quarter of this CCD could be used. CCD #2 and CCD #3 had 
relatively large readout noise. The camera was upgraded in 
2011, at which time these two CCDs were replaced with new 
detectors. CCD # 4 was replaced and swapped with CCD #1. 
New video preamps were added to all four CCDs in order to 
reduce crosstalk and improve noise immunity. 


2.3. Filter 


The SCUSS u filter is similar to the SDSS u band. Figure 2 
displays both SCUSS and SDSS system response functions. 
The SCUSS u response curve includes the filter transmission, 
the CCD quantum efficiency, and the atmospheric extinction at 
a typical airmass of 1.3. The adopted atmospheric extinction is 
the same as that of the SDSS, which is based on the standard 
Palomar monochromatic extinction coefficients but here scaled 
to the elevation of Kitt Peak assuming an exponential scale 
height of the atmosphere of 7000m (Doi et al. 2010). The 
effective wavelength and bandwidth of the filter are respec- 
PrRda PrRda 
[ROA JRA’? 
R(A) is the filter response curve. The effective wavelength and 
bandwidth of the SCUSS u band are about 3538 A and 345 A, 
respectively. The FWHM is about 520A. The effective 
wavelength, bandwidth, and FWHM of the SDSS u filter is 
3562, 385, and 575 A, respectively. The SCUSS u band is 


tively defined as Aeff = and Aeff = where 
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Figure 1. Layout of the CCD array. There are four CCDs: CCD #1, #2, #3, 
and #4. Each CCD has four amplifiers highlighted by dashed lines. Overscan 
regions are also indicated. 


Table 2 
Detector Parameters 


CCD number 4 (in 2 x 2 array) 

CCD gaps 166” in R.A. and 54” in decl. 

CCD size 4096 x 4032 (extra 40 rows of overscan) 
Full well 90,000 electrons /pixel (65,000 ADUs/pixel) 
Amplifier number 4 


~80% at 3500 Å 
7.2 electrons/pixel/hr 


Quantum efficiency 
Average dark current 
Response non-uniformity 2% at u band 
Average readout noise 8.8 electrons 
Readout time ~30 s 

Average gain 1.5 electrons/DN 
Plate Scale 07454 (15 um) 


slightly bluer than the SDSS filter. In the rest of this paper, we 
will use the symbol u* to refer to the SCUSS u band and the 
symbols u, g, r, i, and z for the five SDSS bands. Throughout 
the paper, objects are classified as point-like or extended using 
the SDSS star-galaxy separation unless otherwise specified. 


3. BASIC IMAGE PROCESSING 
3.1. Image Division and Overscan Correction 


90Prime has four CCDs and each CCD is read out by four 
amplifiers. Thus, every exposure file includes 16 FITS 
extensions. We split the file into four smaller FITS images, 
each of which represents one of the four CCDs. Overscan lines 
(40 pixels) are moved to the right side of the image. 

We compute the median of the overscan columns for each 
amplifier and subtract it from the raw frame. There are some 
subtle issues about the overscan that require special care. For 
example, when a bright star is located right beside the overscan 
region, nearby overscan pixels are contaminated. Occasionally, 
some brighter stripes appear within or close to either end of the 
overscan region. The overscan in these regions needs to be 
interpolated or extrapolated. After overscan subtraction, frames 
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Figure 2. Response curves of the SCUSS u* and SDSS u filters. Both curves 


include the atmospheric extinction at an airmass of 1.3 and normalized to their 
maximums. 
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are rotated counter-clockwise by 90° to keep north up and east 
left. These frames are then trimmed to the size of 4096 x 4032 
pixels. 


3.2. Dark Current 


The average dark current of the 90Prime imager is about 7.2e 
per hr. For an exposure of 300s, it is ignorable (about 0.6e) 
relative to the readout noise of about 8.8e, so we did not take 
dark exposures. 


3.3. Bias 


Bias frames are taken before and after the scientific 
exposures. A total of 20 bias frames are obtained each night. 
After combining them, there are large-scale structures, 
especially for CCD #2 and #3 as shown in Figure 3. Counts 
in some peaks of the structures can be more than 10 DN. These 
structures affect the accuracy of photometry significantly, 
especially in u band, because counts of scientific frames in this 
band are low. The bias structures are stable for a long time 
(month-to-month variation is about 0.25 ADU), so we derive a 
median “super” bias by combining all the bias frames taken 
within a month. All overscan-subtracted frames are corrected 
by this “super” bias. 


3.4. Flat-fielding 


Flat-fielding removes the instrumental signatures from raw 
frames using several exposures taken with the telescope facing 
a uniform light source. At the beginning and end of each night, 
a total of 20 dome flats are taken with the dome closed and the 
telescope pointing to a white screen. This screen is illuminated 
by UV lamps of Philips MasterColor Ceramic Metal Halide 
ED-17. Usually, a 6s exposure generates about 20,000 DN on 
the CCD. Although dome flats have extremely high S/Ns, there 
are some disadvantages to applying them in flat-fielding. First, 
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Figure 3. “Super” bias image of four CCDs in 2011 December. The layout of 
these CCDs is the same as that shown in Figure 1. 


lamps are point sources and the scattered light on the screen 
might be uneven in the radial direction. Second, the optical 
path of the light from the screen is not the same as the light 
from night sky. In addition, we also find that the gain of each 
amplifier changes considerably during the whole night. When 
using only dome flats, we find that average sky backgrounds in 
the four amplifiers of each CCD can be quite different. 

Twilight flats are an alternative way and regarded as being 
more uniform than dome flats. They are obtained by observing 
the sky during evening and morning twilight. The brightness of 
the twilight sky changes rapidly and it strongly depends on the 
weather. It is difficult to get enough good high S/N flats. The 
problem of gain variation also persists. 

Raw science images also contain the flat-field characteristics. 
By combining all the science frames with outliers rejected, we 
can obtain a “super” sky flat. This kind of flat has more 
advantages than the two types of flats above. The light of night 
sky “super” flat comes through the same optical path as the 
observed objects. The CCD gain is the average of all science 
frames and it is closer to the real-time variation during the 
night. There are some structures in the “super” flat that present 
a different sensitivity amplitude relative to the dome flats 
possibly due to a different light path. About 150 science images 
on average are observed each night. The average level of the 
sky background in each image is about 150 ADU. As a result, 
the “super” sky flat has an average count of about 23,000 ADU. 
We present the sky flats on 2011 December 28 as an example 
in Figure 4. After applying either the dome flat (higher S/N) or 
super sky flat, we find that there is less than a 0.2% difference 
in the sky background fluctuations. The resulting error due to 
lower S/N of the “super” flat is negligible relative to the large 
sky background fluctuation of about 12.7 ADU and CCD 
readout noise of 8.8 electrons. Thus, we use the “super” sky flat 
to do flat-fielding corrections of the science frames. 
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Figure 4. “Super” sky flat image of four CCDs on 2011 December 28. The 
layout of these CCDs is the same as that shown in Figure 1. 


3.5. Crosstalk 


Crosstalk frequently occurs in multi-channel CCD chip read- 
out. When there are saturated objects in one of the CCD 
amplifiers, it can be best seen as mirror images in other 
amplifiers. The crosstalk usually causes contaminations across 
the output amplifiers at the level of 1:10,000 (Freyhammer 
et al. 2001). The output counts of each amplifier are the sum of 
the true counts from the sky and a small fraction of counts from 
other amplifiers. The crosstalk signal can be either positive or 
negative and should be corrected for high-precision astronom- 
ical photometry. 

The 90Prime camera has four CCDs and each CCD has four 
amplifiers. The crosstalk effect is clearly seen in our SCUSS 
raw images with saturated stars. The crosstalk signal is 
positive. We characterize the effect assuming that it is additive 
and proportional to the number counts of other amplifiers. The 
proportionality coefficients are relatively stable in SCUSS 
images. We correct crosstalk using the following prescription: a 
series of images with bright stars appearing in one of the CCD 
amplifiers is selected; the proportionality coefficients are 
estimated by comparing the mirror signals with their original 
signals; the crosstalk signals in one quadrant are removed with 
the corresponding coefficients of the other three quadrants. The 
average coefficient for the SCUSS images is about 2:10,000. 
Figure 5 illustrates the crosstalk and the performance of its 
correction. The crosstalk signals are clearly seen as mirror 
images (arrows) of a bright star that is marked with a circle. 


4. ASTROMETRIC CALIBRATION 
4.1. Astrometry by UCAC4 


Astrometric solutions are derived by cross-identifying 
objects in science frames with the Fourth US Naval 
Observatory CCD Astrograph Catalog (Zacharias et al. 2013, 
UCAC4). UCAC4 is an all-sky star catalog that is complete in 
R band down to 16 mag. It contains proper motions for most 
stars. An approximate first-order guess of the astrometric 
solution for each CCD is first estimated using UCAC4 objects 
in the same area with positions revised according to their 
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proper motions. Then, the projection center of the telescope can 
be calculated using this rough solution. A more accurate 
astrometric solution is then derived with the projection center 
and a radial second-order correction for the focal plane 
distortion. This distortion term is the same for all frames in 
our survey. 

The left panel in Figure 6 shows the distribution of the 
astrometric errors in the plane of R.A. and decl. differences 
between SCUSS and UCAC4. It includes objects in one 
exposure of a randomly selected field. About 1800 objects are 
matched with UCAC4. The average position offsets are 0” 002 
and 0”004 for R.A and decl., respectively. The lø position 
errors in R.A. and decl. are about 0” 128 and 0” 120. The right 
plot of Figure 6 illustrates the distortion of the focal plane. The 
plate scale varies from the center (0455) to the edge (0” 447) 
of the FOV. The pixel area near the corner is about 3.3% 
smaller than that of the central pixel. 


4.2. External Astrometric Errors 


External astrometric errors are estimated by matching 
SCUSS objects with the UCAC4 catalog. About 165 objects 
on average for each CCD are crossing-identified and the global 
external position error is about 0”13 + 0.02. The mean 
external astrometric offsets and rms errors for R.A. and decl. 
are 


Aa = —0"0012 + 0.0079, A = —0”0015 + 0.0071, 
OAa = 0” 1111 pie 0.0178, ON6 = 0” 1094 pie 0.0182, 


where o corresponds to the 68.3% confidence level. 

Figure 7 displays the external coordinate offset and rms error 
as functions of R.A. and decl. The offsets and rms errors are 
uniform over the entire survey. 


4.3. Internal Astrometric Errors 


Internal astrometric calibration errors are estimated using the 
cross-identifications of UCAC4 sources inside the overlapping 
areas between two adjacent exposures. Each field is observed 
with two exposures dithered by half a CCD size, so every CCD 
frame is covered by four other frames, whose astrometric 
solutions are independently derived. We calculate the internal 
astrometric uncertainties using the objects in common to both 
frames. The global average internal error is about 0” 09 + 0.03. 
The mean internal R.A. and decl. offsets and rms errors over 
the whole survey are 


Aa = 0/0009 + 0.0102, AS = 070004 + 0.0113, 
TAa = 00728 + 0.0295, oas = 070707 + 0.0269. 


5. PHOTOMETRIC CALIBRATION 
5.1. External Calibration 


The SDSS imaging survey has covered more than one third 
of the sky within both northern and southern Galactic caps. It 
provides photometric catalogs of about 5200 square degrees in 
the SGC. More than 21 million objects in this area are 
recorded. The exposure time for each band is about 54 s. The u 
magnitude limit with 95% completeness for point sources is 
about 22.0mag. The photometric calibration accuracy in u 
band is about 1.3%. The ninth data release (Ahn et al. 2012, 
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Figure 5. Left: a SCUSS raw CCD image with visible crosstalk signals. The arrows point to the mirror images of the bright star marked with a circle. Right: the image 


after correcting the crosstalk. 
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Figure 6. Left: distribution of astrometric errors in the plane of R.A. and decl. differences between SCUSS and UCAC4 positions for a randomly selected field. The 
red cross denotes the coordinate origin. Histograms showing the marginal distributions of the R.A. and decl. errors are plotted in the upper and right panels, where the 
red curves show the Gaussian fits. The means and standard deviations of the two distributions are also presented in the plot. Right: pixel scale variation in the detector 
plane. The layout of the CCD array (dashed boxes) is also plotted. Numbers labeling the contours give the pixel scales in arcseconds, which are also indicated by the 


colored surface. 


DR9) is utilized to make photometric calibrations of the 
SCUSS images. 

Because of gain and observational condition variations, 
individual zeropoints must be determined separately for each 
frame. Aperture photometry is performed using DAOPHOT 
(Stetson 1987) on the bias-subtracted and flat-fielded images. 
We choose an aperture radius of about 7/3 (16 pixels), similar 
to that of 7”43 adopted by the SDSS photometric calibration. 
The aperture diameter is about 7 times the typical seeing 
FWHM, which is large enough not to be affected by aperture 
correction. The transformation equation between the instru- 
mental magnitudes uing = —2.5logig ADU to the SDSS 


calibrated magnitude is simply given by 
gt fi 


where c is the instrumental zeropoint, k is the atmospheric 
extinction coefficient, X is the airmass, u — g is the SDSS 
color, and f(u — g) is the color term. The final term of 
f'(u— g, X) is related to the color effect due to the 
atmospheric extinction. The central wavelength of SCUSS u* 
band at an airmass of 1.0, 1.3, and 1.5 is estimated to be 3534, 
3538, and 3542 A. The photometric effect in the color terms for 
most main-sequence stars at different airmasses is less than 


uU = Uins + € + kX + fu 


X%  () 
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Figure 7. SCUSS R.A. and decl. offsets and rms errors relative to UCAC4 as functions of the equatorial longitude and latitude. The two curves in each panel show 


0.3%, estimated using theoretical stellar spectral libraries. So, 
we ignore the second-order color term. Since the standard stars 
are the common SDSS objects within the same area of each 
image, the term solely related to the airmass is constant. Thus, 
Equation (1) can be simplified as 


u* =u — f(u — g) = Uins + C, (2) 


where u is the calibrated SCUSS u-band magnitude and C is 
still termed as the “photometric zeropoint.” The color term is 
estimated as follows: (1) SCUSS instrumental magnitudes are 
calibrated using the photometric magnitudes of common SDSS 
objects without considering the color term; (2) magnitude 
differences between these calibrated SCUSS magnitudes and 
SDSS magnitudes are calculated as a function of the SDSS 
u — g color; (3) an approximate system transformation formula 
is fitted by a two-order polynomial: 


u*— u=- f(u — 8) 
= 0.1459 — 0.1957 (u — 8) 


+ 0.0591 (u — 8}, (3) 


for 0.8 < u — g < 2.7 (see Figure 8). This color term is only 
used for the photometric calibration. Here the SDSS is used for 
the zeropoint only, and the SCUSS photometry on its native 
system is defined to match the SDSS at u — g = 0. We choose 
stars with 16.0 < u < 20.5 to eliminate objects that are 
saturated or have large photometric errors. Equation (3) is 
then applied to transform u to SCUSS u*. Following 
Equation (2), we measure the final zeropoint C by iteratively 
rejecting outliers. 

The CCD gain varies slightly during the observation. 
Furthermore, CCD gains of the four amplifiers did not change 
synchronously. In addition, there is a photometric response 
non-uniformity on the flat-fielded images. It is possibly caused 
by the focal plane distortion and scattered light reflected in the 
optical system (Regnault et al. 2009; Betoule et al. 2013). Thus, 
we perform photometric calibration for each amplifier of each 
CCD. There are about 47 stars in each amplifier on average to 
derive the zeropoint, whose accuracy is estimated to be about 
0.01 mag. The zeropoint here is expressed in units of mag for 
1 ADU s™'. After deriving independent zeropoints for each 
amplifier, we find a small residual pattern in the sensitivity as 


+1o scatter around the average along the x axis. The R.A. is transformed to the range of (—180°, 180°). 


1.0 


05 20 25 3.0 


Figure 8. Magnitude difference between the SCUSS u* and SDSS u as a 
function of the SDSS u — g color. The points are for point-like objects with 
16 < u < 20.5. The error bars show the averages and 1ø errors within different 
color bins. The red curve is the fitted transformation equation as shown in 
Equation (3). 


Figure 9. Non-uniformity pattern used to correct the flat-fielded and calibrated 
images. The layout of these CCDs are the same as in Figure 1. 
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derived from objects observed in different positions in the field 
as shown in Figure 9. We use this residual map to further refine 
the flat-field. 
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Figure 10. Left: average zeropoints of the four amplifiers for each CCD as a function of time. The data points correspond to the CCD images, which are sorted by the 
observation time. The zeropoints are in units of mag for 1 ADU s™'. The data for CCD #1, #2, and #3 are shifted by 2.0. The dashed lines separate the years of 
observations. Right: relative zeropoints compared with the above average for each amplifier of CCD #1 as a function of time. The data for three amplifiers, Amp. #1, 


#2, and #3, are shifted by 0.2. 


5.2. Photometric Response Differences of the Four Amplifiers 


The average zeropoints of four amplifiers are shown in the 
left panel of Figure 10 as a function of time. Usually during the 
night, we scan the sky from east to zenith and then to west. 
The airmass and corresponding atmospheric extinction first 
decrease and then increase. Therefore, zeropoints present 
nightly variations. From this figure, it can be seen that the 
weather conditions over 2010 are worst. The right panel of 
Figure 10 shows the relative zeropoint variation with time for 
each amplifier of CCD #1 compared with the average 
zeropoint of all amplifiers. The relative zeropoint variations 
do not keep the same value and sometimes two amplifiers 
present opposite variations. The detectors have bad qualities in 
2010, so the relative zeropoint changes more notably with time. 
Thus, it is critical that we obtain the photometric solutions 
independently for the four amplifiers. 

Table 3 gives the zeropoint scatter of four amplifiers for each 
year and each CCD. The zeropoint differences among 
amplifiers originate from gain variation, photometric response 
non-uniformity as mentioned before, and the photometric 
calibration error. The overall zeropoint scatter of four 
amplifiers around the averages are 0.012 + 0.006, 0.016 + 
0.007, 0.015 + 0.009, and 0.011 + 0.006 for CCD #1, #2, 
#3, and #4, respectively. The general response differences of 
four amplifiers in all CCDs are less than 1.5% except CCD #3 
of 2010, which is about 2.8%. 


5.3. Internal Calibration 


Each SCUSS field is observed twice and overlaps with the 
adjacent fields. Thus, each CCD image can be calibrated using 
the common objects cross-identified in surrounding images. 
The calibrating process begins with balancing the zeropoints of 
images covered by the SDSS and then transferring the 


Table 3 
Average Zeropoint Scatter of Four Amplifiers for Each CCD 

coD #1 #2 #3 #4 

year 

2010 0.011 + 0.007 0.015 + 0.009 0.028 + 0.009 e 

2011 0.013 + 0.006 0.015 + 0.006 0.011 + 0.006 0.010 + 0.006 
2012 0.013 + 0.006 0.013 + 0.006 0.010 + 0.006 0.011 + 0.006 
2013 0.010 + 0.006 0.020 + 0.007 0.017 + 0.006 0.011 + 0.006 
All 0.012 + 0.006 0.016 + 0.007 0.015 + 0.009 0.011 + 0.006 


photometric solutions to the images out of the SDSS footprints. 
The calibration iterates until the whole grid of photometric 
solutions finally converges. Figure 11 shows the magnitude 
difference distribution between measurements for objects that 
are observed twice. The histograms in black and red are the 
distributions of the same objects with the SCUSS u-band 
magnitude calibrated externally and internally, respectively. 
The dispersion by external calibrations is about 0.028 mag, 
while that by internal calibrations is about 0.025. It seems the 
internal calibration is at least as good as or even a little better 
than the external calibration. 


6. IMAGE QUALITY STATISTICS 


Both weather and instrumental status affect the image 
quality. The camera was regularly updated after each observa- 
tion season, so its performance was improved gradually. The 
weather and night sky conditions strongly affect the imaging 
depth of the survey. We can measure the characteristics tracing 
the image quality, such as airmass, sky background brightness, 
seeing, and photometric zeropoints. 
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objects observed more than twice. The black histogram is externally calibrated 
by the SDSS catalog, while the red one is internally calibrated. Only bright 
objects with 16 < u < 19 are selected. 


6.1. Seeing 


The seeing is estimated by the FWHM measurements of 
isolated and bright point sources. The seeing distribution is 
presented in Figure 12(a). The best seeing is 1”2 and the 
overall median seeing is about 2”0. The Bok telescope is 
located in the trough between two peaks of the mountain, so the 
wind speed is usually larger than other places on the same site, 
which has an effect on the seeing. 


6.2. Airmass 


Figure 12(b) shows the airmass distribution of all survey 
images. The median airmass is about 1.28, and this is used to 
determine the typical u-band filter response as shown in 
Section 2.3. 


6.3. Photometric Zeropoint 


The variation of the photometric zeropoint mainly reflects 
the change in atmospheric transparency. The distribution of the 
photometric zeropoints is plotted in Figure 12(c). The mean 
photometric zeropoint is 23.81mag for 1ADUs~'. The 
median zeropoints for airmasses of 1.0, 1.2, and 1.4 are 
23.93, 23.85, and 23.75, respectively. There are some images 
with low zeropoints, most of which were observed when the 
weather was cloudy. 


6.4. Sky Brightness 


The sky background brightness is an important parameter to 
quantify a ground-based observation station. Compared with 
other light sources, the artificial light pollution from nearby 
cities is more serious at Kitt Peak. The distribution of the 
u-band sky brightness at Kitt peak is shown in Figure 12(d). 
There are some images taken when the moon is above the 
horizon. The average night sky brightness of all observations 
is about 22.05 mag arcsec. The moonless median sky 
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background at zenith is about 22.37 mag arcsec *, which 


is comparable to that of the Apache Point site 
(22.1 mag arcsec *). Note that the calibrated sky brightness 
uses the frame zeropoints that implicitly include atmospheric 
extinction, so the actual sky brightness is darker because the 
u-band atmospheric extinction coefficient is about 0.5, as 
estimated with the observations taken during photometric 
nights. 


7. IMAGE STACKING 
7.1. Quality Control 


Most SCUSS fields are observed under good weather and 
moderate seeing conditions. There are some cases with bad 
image quality: (1) high sky background due to observations 
during astronomical twilight or when the moon was up; (2) 
large seeing due to strong wind; (3) low atmospheric 
transparency due to cirrus in the FOV; (4) bad focus of the 
CCD camera. In 2013, we spent one and half observation runs 
re-observing most bad-quality fields. 

To ensure homogeneity of the imaging depth and complete- 
ness of the SCUSS survey, we keep only the images with 
seeing <3”0, sky background ADU <500, and photometric 
zeropoint >22.56 mag. About 92.6% of the survey area is 
covered by these images. For the remaining area, we take 
images with seeing <3”0 (3.6% area), regardless of the sky 
brightness and photometric zeropoint. If none are available, we 
take all remaining images, which covers about 3.8% of the 
total area. 


7.2. Image Resampling and Stacking 


Based on the central coordinates of each field, we stack 
related single-epoch images to form a combined image. We 
first project and resample the single-epoch images to a grid 
with a fixed pixel scale of 0”454. The grid has 8640 x 8200 
pixels, covering a sky area of 12090 x 12034. If a single- 
epoch images contributes less than less than 128 x 128, it is 
not included in the stacking process. For each pixel of the grid, 
there are more than four related pixels in each single-epoch 
image. These pixels are regarded as having the same size after 
being flat-fielded. We calculate the fractions of these pixel 
areas that are covered by the grid pixel and sum them to 
conserve flux. 

We subtract the sky backgrounds from the resampled 
images. Their photometric zeropoints are converted to linear 
flux weights. The remaining signal after removing backgrounds 
is weighted by these weights and then co-added. If there are 
more than three pixels involved in stacking, a cosmic-ray 
rejection is implemented using a sigma clipping algorithm. We 
redo the flux calibration for each stacked image with the SDSS 
DR9 catalog to derive the final photometric zeropoint. This 
zeropoint is approximately 29 mag for 1 DN. In addition, a 
mask image for each field is also generated. Each mask pixel 
presents the number of images that are actually involved in the 
flux co-adding. The mask value is reduced by one for each 
epoch in which a pixel is bad, saturated, or blank. Figure 13 
shows a typical stacked image and its mask image. Most of the 
stacked image has two exposures, an overlapping area has more 
than three exposures, and there are some small holes located in 
the CCD gaps that are blank. Some of the CCD gaps have only 
one exposure, so the depth is about 0.75 mag shallower in those 
locations. 
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Figure 12. Histograms of seeing (in arcsec), airmass, photometric zeropoint (in mag for 1 ADU s_'), and sky background (in mag arcsec *) for all SCUSS u-band 


observations. 


Figure 13. Examples of a stacked image (left) and its mask image (right). The central coordinate is a = 44°7536, 6 = 6° 1074. North is up and east is left. A galactic 
cluster of Abell 400 (marked as a plus symbol in the left plot) happens to appear in this area. The mask image shows the exposure numbers. There are two black boxes 
presenting no observations. Most of the area has two exposures, a small part shows more than two exposures, and some CCD-gap areas have only one exposure. The 
small horizontal lines are caused by the missing CCD rows that were not recorded by the camera controller during 2011 and 2012. 
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Table 4 
Aperture Radii for the Aperture Photometry 


No. Radius in Pixels Radius in Arcseconds 
1 3 1.36 
2 4 1.82 
3 5 2.27 
4 6 2.72 
5 8 3.63 
6 10 4.54 
7 13 5.90 
8 16 7.26 
9 20 9.08 
10 25 11.35 
11 30 13.62 
12 40 18.17 


8. PHOTOMETRY 


The SCUSS photometric pipeline generates comprehensive 
catalogs using multiple photometric techniques so that different 
users can choose the one that suits their needs. Aperture, 
automatic aperture, point-spread function (PSF), and model 
photometry are applied to both SCUSS stacked and single- 
epoch images. The model photometry is consistent with the 
SDSS model photometry that utilizes the SDSS r-band model 
shape parameters to measure brightnesses on the SCUSS 
u-band images. 


8.1. Source Detection 


Source detections for astronomical science images are never 
complete, especially at the fainter magnitude end. Morpholo- 
gies of extended sources in u band are more fragmented and 
diffuse than those in other optical bands. Therefore, many 
fainter sources with low surface brightness might be missing. 
In addition, most of the science projects based on SCUSS data 
also need to use deeper and redder other bands from other 
large-scale surveys. More than 3/4 of the SCUSS area is 
covered by the SDSS. Our sources include both SDSS detected 
objects and detections unique to our SCUSS images. For the 
area not covered by the SDSS, we perform photometry for only 
SCUSS detections. The resulting catalog within these areas will 
be useful to match with other wide imaging surveys (e.g., Pan- 
STARRS). 

SDSS objects with any one of the PSF, Petrosian, model, and 
CModel ugriz-band magnitudes brighter than 23.5 mag are 
selected. SExtractor is used to detect objects in SCUSS stacked 
images, and these are matched with the SDSS objects. The 
mismatched objects are rematched with SDSS catalogs with 
magnitudes fainter than 23.5 mag in order to find missing 
SDSS objects. The rest of the objects are SCUSS unique 
detections. 


8.2. Aperture Photometry by DAOPHOT 


Circular aperture photometry is a simple procedure to 
measure magnitudes of sources. The core code of the aperture 
photometry in DAOPHOT is utilized to measure aperture 
magnitudes (Stetson 1987). We use 12 apertures with radii 
ranging from 14 to 18” (from 3 to 40 pixels; see Table 4). 

It is important to apply aperture corrections to aperture 
magnitudes due to flux loss within a finite aperture size. For 
smaller apertures and larger seeing, this is especially important. 
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Figure 14 presents growth curves under different seeing 
conditions. The growth curve is calculated as the magnitude 
difference between 1 of 12 apertures and the 73 radius 
aperture as function of aperture radius. This reference aperture 
is the one used for photometric flux calibration. We choose 
isolated and point-like objects with photometric errors in all 
apertures less than 0.05 mag to derive the growth curve. Outlier 
objects are eliminated when the median growth curve is 
calculated. 


8.3. Automatic Photometry by SExtractor 


SExtractor (Bertin & Arnouts 1996) provides precise 
magnitudes of sources using automatic aperture photometry, 
generating a so-called “automatic magnitude,’ which is 
motivated by the Kron algorithm (Kron 1980).° An elliptical 
aperture is automatically determined for each object to integrate 
the flux. Although most of the flux is expected to lie within the 
elliptical aperture and the flux loss should be almost 
independent of the source magnitude, we discover that in fact, 
the flux loss does change with both the source magnitude and 
seeing. The magnitudes of objects that are brighter or objects 
observed with worse seeing have more corrections. We 
consider an equivalent circle whose area is equal to that of 
the ellipse. Thus, the radius of this circular aperture is equal to 
the root of the product of the semi-major and semi-minor axis 
lengths. We call this circular radius the characteristic radius, 
which has the same meaning as the aperture size in normal 
aperture photometry discussed above. Figure 15 presents the 
aperture corrections for the automatic magnitudes with 
different seeing. The black dots are point-like and isolated 
objects detected by SExtractor. They lie along with the growth 
curve of the aperture correction described in the previous 
section. Thus, an interpolation from the curve is good enough 
for correcting automatic magnitudes. 

We also perform Petrosian-like photometry on SCUSS 
images by SExtractor. The Petrosian aperture is very similar to 
the Kron one. They share the same position angle and 
ellipticity. The Petrosian aperture radius is determined by the 
ratio of the isophotal brightness at a certain radius and average 
surface brightness within this radius. The ratio is set to be 0.2 
and the corresponding Petrosian radius is larger than the Kron 
radius. The aperture corrections for Petrosian magnitudes can 
be estimated the same way as mentioned above. 


8.4. PSF Photometry 


The PSF is obtained using PSFEx’ (Bertin 2011). The form 
of the PSF in PSFEx is expressed as a linear combination of 
basis vectors. The pixel basis is selected in our PSF modeling. 
The spatial PSF variation on the focal plane usually shows a 
smooth profile, which can be modeled by a low-order 
polynomial. For SCUSS single exposures, a second-order 
polynomial is good enough to describe the PSF variation over 
the CCD plane. For SCUSS stacked images, a seventh-order 
polynomial is used to describe the complexity of the image 
quality. 

A code is specially designed to measure the PSF magnitudes 
at known object positions. The code takes the SCUSS image, 
SDSS and SCUSS-only object positions, the bad-pixel list, and 


6 https: //www.astromatic.net /software /sextractor 
7 http: //www.astromatic.net/software/psfex 
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Figure 14. Growth curves for aperture corrections at different seeing. The black pluses are growth curves for about one hundred isolated and point-like objects. The 
red lines with error bars are the medians and standard deviations of the corrections for different apertures. The last panel shows the median growth curves in the 


previous three panels in order to present the impact of seeing. 


0.1 FS 
—O.0F 3 


seeing: 


seeing: 2.0 


10 12 14 4 
r (pb 


Figure 15. Difference between the automatic and circular aperture magnitudes of 7”3 radius as a function of the characteristic radius. Black dots are isolated point 
sources. Growth curves for corresponding circular aperture corrections shown in Figure 14 are overlaid in red. The cyan curves are the cubic spline interpolation of 


those growth curves. 


the PSF model derived by PSFEx as inputs and then outputs the 
Gaussian-fitted position, local sky background and its error, 
local PSF FWHM, fitted PSF integrated flux and its estimated 
error, and a flag tagging the status of each object. Figure 16 
shows the flowchart of the PSF photometry code. The local sky 
background for each object is measured using the pixels at 
r> 71.5 x FWHM. Here, the FWHM is the FWHM of the 
local PSF model. Outliers, such as cosmic rays and signals 
from real objects, are iteratively rejected by a sigma-clipped 
algorithm. After the sky background is subtracted, the object 
position on the CCD is fitted by a two-dimensional Gaussian 
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function. The position is allowed to shift because the 
coordinates in SDSS redder bands are slightly different from 
those in the SCUSS u band due to atmospheric refraction and 
star proper motions. The pixels centered at the initial position 
with r < 2.5 x FWHM are considered in the calculation. If the 
new fitted position is more than 2.0 pixels away from the old 
one, the following PSF photometry will be performed at the 
original position. The PSF model is interpolated to the same 
pixel scale of the CCD image and fitted to the fluxes of the 
object pixels within r= 2.0 x FWHM. Flags are also 
provided by our PSF photometry to show the reliability. They 
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Figure 16. Flow chart of our PSF photometry pipeline. First, the PSF model is extracted by combining SExtractor and PSFEx. Second, the coordinate for each object 
is fitted by a Gaussian function. Third, the image of each object with the sky background subtracted is fitted by the local PSF profile. Finally, the photometric results 


together with the fitted positions and flags are given as outputs. 
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Figure 17. Comparisons of our PSF photometry with DAOPHOT. Left is the comparison for common bright objects with 16 < u < 18 and right is the one for 


common faint objects with 21 < u < 22. 


are coded in decimal and expressed as a sum of powers of 2: (1) widely used package for accurate stellar photometry designed 
CCD artifacts; (2) bad pixels; (4) including saturated pixels; (8) to deal with crowded fields. We use DAOPHOT to find objects 
contaminated by neighbors; (16) near image edges. in one single-epoch image and at the same time give the PSF 

We compare our PSF measurements with one of the popular magnitude measurements. Then, our code performs PSF 


PSF photometry software packages (DAOPHOT), which is a photometry for those objects with the PSF model derived by 
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PSFEx. The magnitude comparison is shown in Figure 17. Two 
groups of objects are chosen: a brighter one with 16 < u < 18 
mag (left in Figure 17) and a fainter one with 
21 < u < 22 mag (right in Figure 17). The scatter of the 
PSF magnitude differences between our code and DAOPHOT 
are 0.009 mag for the bright group and 0.021 for the faint 
group, respectively. Our PSF photometry is consistent with 
DAOPHOT, although the PSFs used by these two techniques 
are derived in different ways. 


8.5. Model Photometry 


The SDSS also provides a type of model magnitude, or 
“modelMag,” which can measure the unbiased colors of 
galaxies in the absence of color gradients through equivalent 
apertures in different bands. The “modelMag” is generated by 
choosing a shape model, either a deVaucouleurs or exponential 
profile, based on its best-fit likelihood in the SDSS r band, then 
convolving them with seeing in other bands, and finally forcing 
magnitude measurements with the same aperture shape. 

We divided SDSS objects into two groups, point sources and 
extended sources, based on the SDSS classification. For 
extended sources, we construct their theoretical 2D models 
with the effective radii, axis ratios, and position angles from the 
SDSS measurements. These models are convolved with local 
PSF profile derived by PSFEx. For point sources, we use local 
PSF from PSFEx directly as their models. In addition, the 
extended sources with small sizes or low brightnesses (effective 
radii less than 0.5 arcsec or SDSS r-band magnitude fainter 
than 23.5 mag) are also treated as point sources. 

After the models are constructed, model amplitudes are 
calculated as the ratios of the models and raw SCUSS fluxes. 
The pixels within r = 1.0 x FWHM are considered in the 
calculation, which is optimized for faint SDSS extended 
sources. According to the amplitudes, we then compute the 
corresponding deVaucouleurs and exponential magnitudes. 
Therefore, the “modelMag” in SCUSS u band can be estimated 
by the SDSS r-band deVaucouleurs and exponential profiles. 
The SCUSS model magnitude is aperture-corrected to make the 
magnitudes of bright point sources (16 < u < 20) equal to the 
SDSS u-band model magnitudes. 


8.6. Co-added Photometry 


The stacked images are composed of single-epoch images 
taken with different observational conditions. The seeing varies 
between each single-epoch image, which makes the PSF profile 
of the stacked image quite complicated. Subsequently, the 
fraction of flux loss due to finite apertures for objects in 
different parts of the stack images might be different. In 
addition, the PSF profile cannot be perfectly determined unless 
the seeing of related single-epoch images is similar enough. On 
the other hand, the PSF profile of a single image varies 
smoothly and it is much easier to model. We can first perform 
photometry for an object observed with different exposures and 
then co-add its fluxes to generate different co-added magni- 
tudes at the catalog level. 

The aperture photometry for single-epoch images is obtained 
by DAOPHOT with 12 apertures as defined before. Aperture 
corrections are applied to the resulting aperture magnitudes. 
Corresponding corrected fluxes are weighted by the errors and 
averaged to generate the co-added aperture magnitudes. 
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The PSF photometry for single-epoch images is performed 
using our PSF fitting code. The PSF magnitudes are also 
aperture-corrected by comparing with the aperture magnitudes 
in a 7”3 radius. The co-added PSF magnitudes are calculated 
by averaging the fluxes weighted by their errors. 

We perform the model photometry for single-epoch images 
and measure the exponential and deVaucouleurs magnitudes. 
We adopt the model magnitude with a higher SDSS r-band 
likelihood. Aperture corrections are applied to make model 
magnitudes equal to the SDSS model magnitudes in the case of 
unresolved objects. The co-added model magnitude is 
calculated by averaging the model fluxes weighted by their 
errors. 


9. SOME PHOTOMETRIC COMPARISONS 


The photometric methods described above are applied to 
stacked and single-epoch images. Figure 18 illustrates the 
general photometric scheme. The detections are based on both 
SDSS and SCUSS images. We perform aperture, PSF, model, 
and automatic photometry for the stacked images. Note that the 
automatic magnitudes are only based on objects solely detected 
by SExtractor. We also obtain co-added aperture, PSF, and 
model magnitudes from photometry of single-epoch images. 
Since point sources can be best modeled while extended 
sources are much more complicated, the photometry of point- 
like objects is better suited for comparisons of our different 
photometric methods. Thus, the comparisons of point sources 
are presented for most cases below and, for comparison with 
SDSS, we correct the magnitude with the SCUSS/SDSS 
color term. 


9.1. Photometry of Stacked Images 


We compare the different SCUSS photometric measure- 
ments for stacked images with the SDSS PSF or model 
magnitudes of point sources in Figure 19. The SCUSS model 
magnitude is compared with the SDSS u-band model 
magnitude, while other magnitudes are compared with the 
SDSS PSF magnitude. We choose the 5 pixel (227) to present 
the aperture magnitude in this figure. For the smallest apertures, 
the aperture magnitudes might be problematic due to 
inhomogeneous image quality in some stacked images. Table 5 
gives the magnitude offset and scatter in two magnitude 
intervals: 16 < u < 19 mag and 20.5 < u < 21.5 mag (i.e., 
~21 mag). The automatic magnitude is best among all 
photometric magnitudes for stacked images, with a scatter of 
0.033 for bright stars and 0.178 for faint stars at u ~ 21.0 mag. 


9.2. Co-added Photometry for Single Images 


The comparisons of the co-added PSF, aperture, and model 
magnitudes with the SDSS PSF or model magnitudes are 
shown in the right panels of Figure 19 and the last three rows of 
Table 5. The co-added PSF magnitude performs the best for 
point sources among all the magnitude types. The scatter is 
0.033 for brighter sources and 0.174 for fainter sources at 
u ~ 21.0 mag. The co-added aperture magnitude is also 
adequate if a proper aperture radius is considered. The aperture 
size should be chosen according to the object type, object 
brightness, signal-to-noise requirement, etc. 


THE ASTRONOMICAL JOURNAL, 150:104 (17pp), 2015 October 


Stacked images 


Single images 


ZOU ET AL. 


SDSS DR10 
coordinates 


Source detection: SCUSS only + SDSS objects 


PSF 
photometry 


Automatic 
photometry 


Magnitudes 
for stacked images 


Model 
photometry 


Aperture 
photometry 


Co-added magnitudes 
for single-epoch images 


Figure 18. SCUSS photometric scheme. The output magnitudes include the automatic, Petrosian, aperture, and PSF magnitudes for the stacked images and co-added 


aperture and PSF magnitudes for the single-epoch images. 


9.3. Comparisons with the CFHTLS Deep u Band 


The Canada—France-Hawaii Telescope Legacy Survey 
(CFHTLS; Astier et al. 2006) used the wide-field optical 
imaging camera on CFHT to obtain deep multicolor photo- 
metry over wide areas over 5 years. The photometric system is 
similar to that of the SDSS, except the u filter. As one 
component of the CFHTLS, the wide synoptic survey covers 
155 square degrees in four patches down to i= 24.5. The u- 
band magnitude limit in the wide survey reaches 25.2 mag, 
which is much deeper than the SCUSS u band. Thus, the 
catalogs of the wide survey can be regarded as a reference 
when comparing the SCUSS photometry with that of the 
SDSS. There are two CFHTLS wide regions, W1 and W4, 
overlapping with the SCUSS footprints. The W4 field 
(œ = 22613™185, 6 = +01419™00%) has an area of 25 deg” 
fully covered by the SCUSS and it is located at a higher 
declination, where the SCUSS data quality is more typical than 
it is in the W1 field. We make some photometric comparisons 
using the W4 catalogs. 

Figure 20 shows the photometric comparisons of point 
sources between the SCUSS and SDSS with the CFHTLS wide 
data as a reference. The left panels in this figure present the 
magnitude difference between the SDSS (in blue points) or 
SCUSS u-band magnitude (in red points) and the CFHTLS 
automatic magnitude as a function of the CFHTLS u magnitude. 
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The SDSS magnitude is the PSF magnitude and the SCUSS 
magnitudes are automatic, co-added PSF, and co-added aperture 
(227) magnitudes from top to bottom, respectively. These three 
types of SCUSS magnitudes are considered to be the best flux 
measurements for point sources. The number of objects for the 
automatic magnitude is less than those of the other two 
magnitude types (mainly at the faint end) due to different source 
detection by SExtractor itself. Both the SCUSS and SDSS 
u-band magnitudes are converted to the CFHTLS photometric 
system by the color term as indicated in the CFHTLS webpage’: 
UCFHT = Uspss/scuss — 9.241 (Uspss/scuss — &spss): The 
SDSS scatter for brighter sources is similar to the SCUSS 
scatter, while it is much larger than the SCUSS scatter at the 
faint magnitude ends. The histograms in the right panels of 
Figure 20 show the magnitude difference distributions with 
22 < uceyr < 23. The scatter for the SCUSS and SDSS is also 
presented in these panels. The photometric accuracy of the 
SCUSS is much better than that of the SDSS for fainter sources 
due to deeper imaging. 


9.4. Photometry for Extended Sources 


Extended sources are much more complicated, especially in 
ultraviolet bands, since their morphologies look much more 


8 http: //cfht.hawaii.edu /Instruments /Imaging /MegaPrime/ 


generalinformation.html 
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Figure 19. Magnitude differences between different SCUSS magnitudes and the SDSS PSF or model magnitudes of points sources as functions of SCUSS 
magnitudes. The first-row panels and the rest of the left panels are comparisons of automatic, Petrosian, PSF, and aperture magnitudes for the stacked images. The 
right panels from the second row to bottom are comparisons of co-added PSF, aperture, and model magnitudes for the single-epoch images. The aperture magnitude 
shown here is the one with an aperture radius of 4”5. The two lines above and below zero show +1ø scatter around the average. 


Table 5 
Magnitude Offset and Scatter of Point Sources for the SCUSS Magnitudes 
Compared with the SDSS PSF or Model Magnitudes 


Magnitude 16<u< 18 20.5 <u < 21.5 

Offset Scatter Offset Scatter 
Automatic (stacked) 0.001 0.033 0.008 0.178 
Petrosian (stacked) 0.002 0.033 0.065 0.200 
PSF (stacked) —0.003 0.047 0.025 0.175 
Aperture (stacked) —0.002 0.044 0.017 0.176 
Model (stacked) —0.006 0.059 0.004 0.177 
Co-added PSF —0.004 0.033 0.021 0.174 
Co-added aperture 0.003 0.034 0.010 0.177 
Co-added model —0.003 0.049 0.003 0.175 


fragmented than in redder bands. By comparing different 
magnitude measurements on stacked images with the SDSS 
model magnitude, we calculate a scatter at u ~ 21 of about 0.2, 
0.31, 0.4, 0.23 for automatic, aperture (2"27), PSF, and model 
magnitudes, respectively. The co-added photometric magni- 
tudes give similar results. The automatic magnitudes seem to 
be the best for extended sources, but they only describe the 
brightnesses of SCUSS-detected objects. The model and 
aperture magnitudes are also adequate. The PSF photometry 
fails to measure the magnitudes of galaxies. 


9.5. Guidelines to use Magnitudes 


The choice of photometric technique is dependent on the 
science one wants to do. Here we present some general 
guidelines. More than 90% of the stacked images are 
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assembled from single-epoch images with consistent image 
quality. The photometry on these images is as good as the 
corresponding co-added photometry. However, as the back- 
ground noise in single-epoch images is larger than that in 
stacked images, there are about 23% more objects with 
available magnitude measurements in stacked images than 
with available co-added magnitudes. 

The automatic magnitude performs as well as the PSF and 
model magnitudes because it can adaptively fit elliptical 
apertures to both point and extended sources with similar flux 
loss. It can be regarded as a universal magnitude for both point- 
like and extended sources. But unlike other photometric 
methods, the automatic photometry is based on the objects 
detected on SCUSS images, which are about 35% of SDSS 
objects with available SCUSS u-band fixed-parameter mea- 
surements. For point sources like quasars and stars, the PSF 
and aperture magnitudes with appropriate aperture sizes are 
recommended. The aperture size needs to be determined based 
on scientific objectives. For nearby galaxies with extended 
morphological structures, the automatic magnitude and aper- 
ture magnitude are good choices. The SCUSS model 
magnitude is defined the same as the SDSS “modelMag.” 
When combining magnitudes of other SDSS bands to measure 
the colors of extended sources, it is better to use the model 
magnitude. 


10. SUMMARY 


The SCUSS survey is a wide-field u-band sky survey in the 
Southern Galactic cap. The survey used the Bok telescope on 
Kitt Peak and the filter is close to the SDSS u band. The survey 
observations were completed by the end of 2013 and the total 
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Figure 20. Photometric comparisons of point sources between the SCUSS (in red) and SDSS (in blue) with the CFHTLS W4 data as a reference. The SDSS magnitude 
is the PSF magnitude and the CFHTLS magnitude is the automatic magnitude. The left panels are the magnitude differences between the SCUSS or SDSS and the 
CFHTLS as a function of the CFHTLS u-band magnitude, and the right ones are the normalized distributions of the differences at the faint magnitude end (the hatch 
area: 22 < u < 23 mag). From top to bottom, they are comparisons of three types of SCUSS magnitudes: the automatic, co-added PSF, and co-added aperture 
magnitudes, respectively. Only about 20,000 objects are randomly selected in this figure in order to avoid crowding. The colored text in the right panels gives the 
magnitude scatter of the SCUSS and SDSS, which corresponds to the 68.3% confidence level. 


area is about 5000 deg”. This paper describes the detailed data 
reduction dedicated to the survey, including basic image 
processing which has some special features related to the 
detectors, astrometric and photometric calibrations, and photo- 
metry. The general astrometric error is about 0713. The 
SCUSS photometric calibration is tied to the SDSS catalogs 
and is performed for each amplifier of a CCD due to gain and 
weather variations. 

We apply different photometric techniques to the stacked 
images, including automatic photometry by SExtractor, 
aperture photometry by DAOPHOT, PSF photometry, and 
model photometry. Our PSF photometry with more control- 
lable parameters is consistent with DAOPHOT. The model 
photometry is similar to the SDSS “ModelMag,” which uses 
the SDSS r-band model-derived shape parameters and SCUSS 
PSF profiles to make consistent and unbiased model magnitude 
measurements. We perform photometry on stacked images and 
also on single-epoch images, from which co-added photometry 
is derived. More than 90% of stacked images are assembled 
from single-epoch images with consistent quality. There are 
about 23% more objects with available magnitudes on stacked 
images than with available co-added magnitudes. The 
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photometry on these stacked images is as good as the co- 
added photometry. However, for the rest of the stacked images, 
their photometry is worse than the co-added photometry due to 
uneven image quality. 
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